Geometric nature of the environment-induced Berry phase and geometric dephasing 
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We investigate the geometric phase or Berry phase (BP) acquired by a spin-half which is both 
subject to a slowly varying magnetic field and weakly-coupled to a dissipative environment (either 
quantum or classical). We study how this phase is modified by the environment and find that 
the modification is of a geometric nature. While the original BP (for an isolated system) is the 
flux of a monopole-field through the loop traversed by the magnetic field, the environment-induced 
modification of the BP is the flux of a quadrupole-like field. We find that the environment-induced 
phase is complex, and its imaginary part is a geometric contribution to dephasing. Its sign depends 
on the direction of the loop. Unlike the BP, this geometric dephasing is gauge invariant for open 
paths of the magnetic field. 

PACS numbers: 03.65.Vf, 03.65.Yz, 85.25.Cp 
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Introduction. The Berry phase (BP) is a fundamen- 
tal quantum-mechanical phenomenon related to the adi- 
abatic theorem. Berry Q showed that the phase acquired 
by an eigenstate of a slowly varying Hamiltonian Hit) is 
related to the geometric properties of the loop traversed 
by H{t). In the presence of dissipation the condition of 
adiabaticity and the existence of the Berry phase require 
careful analysis. The widespread criterion of adiabaticity, 
based on a comparison of the rate of change of the Hamil- 
tonian with the gap in the spectrum should be modified 
to involve the matrix elements of the system-environment 
coupling, cf . .2;] . Here we study the interplay of the vary- 
ing field and the dissipation, analyzing the BP in the 
limiting case of weak system-environment coupling. This 
analysis is relevant to the recent and proposed experi- 
ments to manipulate quantum two- level systems (qubits). 
Our findings could be tested in solid-state qubits, such 
as superconducting nanocircuits 0, 0, 0, • 

Berry yij considers a two-level spin-half system in a 
magnetic field Q, which is varied slowly along a closed 
path: H sp i n — — iB(t)<r. The rate of the field's change 
is characterized by the time to complete the loop, £p. In 
the adiabatic limit, Btp 3> 1, the relative phase acquired 
by the eigenstates is a sum <f> — § \B(t)\dt + $bp of the 
dynamical and Berry phases. The latter is geometric, 
it depends on the geometry of the loop but not on the 
details of its traversal (for an isolated spin-half it is given 
by the solid angle subtended by the loop B(i)). In the 
spin language, the evolution is a rotation of the spin by 
an angle $ about B. 

If the spin is not isolated, the dynamics are more com- 
plicated. For a static field, B, dissipation induces en- 
ergy and phase relaxation processes (with the time-scales 
Ti, T2 respectively) and a Lamb- like shift of the level 
splitting, <5i?Lamb, which modifies the dynamical phase. 



(a) field 5a(B) 

(solid arrows) 

\B 7 



(b) quadrupole-like field 
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FIG. 1: (a) The field 5a(B), see eq. (£3, in th e vicinity 
of an open path of B(t) (the dotted line). The environment- 
induced modification of the Berry phase, <5"1>bp, is the integral 
of 5a(B) along this path. For such open-paths, the geomet- 
ric dephasing is gauge-invariant (unlike the phase), (b) The 
quadrupole-like field Sh(B), see eqs. I14I15L (Sh(B) is cylin- 
drically symmetric about the B 2 -axis). For a closed path of 
B(t), <5$bp is given by the flux of 5h(B) through the path. 



Strong dephasing masks the BP, however it can be ob- 
served if B ;§> T 2 _1 2]. In this weak-coupling limit one 
can carry out a BP experiment slowly enough for non- 
adiabatic effects to be ignored (tp 3> B^ 1 ) while ensur- 
ing it is fast enough that dephasing has not destroyed all 
phase information (tp < T2). 

In this letter we show that the coupling to the en- 
vironment modifies the BP when the magnetic field is 
slowly varied. In the adiabatic limit this modification, 
<5$bp, is geometric (see Fig. and complex Its real 
part is an environment-induced BP, thus the total BP is 



$bp + Re[(5$ B p] (where $£ p is the BP for an isolated 
spin-half). Its imaginary part, Im[5$Bp]; is a geomet- 
ric correction to dephasing, whose sign depends on the 
direction of the loop. The magnitude of coherency will 
thus depend on the sign of the loop's winding number. 
Notably Im[i5 c f ) Bp] is gauge-invariant for open as well as 
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closed paths of the i?-field. 

Specifically we study the environment-induced BP for 
an arbitrary path B(i), allowing us to analyze its geo- 
metric nature. For an isolated spin-half the Berry phase 
is given by the flux through the closed loop B(f) of the 
field bjn(B) of a charge-one monopole at the origin in B- 
space yj: we show that the environment-induced modi- 
fication may be interpreted similarly and find the corre- 
sponding field distribution <5b(B) in the B-space. This 
field (and thus <5<&bp) scales quadratically with the spin- 
environment coupling, similarly to the dissipative rates, 



rp-l rp- 



If the environment-induced field fluctuates 



along a single direction, z, the field <5b(B) is axially sym- 
metric, has no axial component, 5b v = 0, and the angular 
distribution of a quadrupole. The field <5b(B) is given by 
a sum over frequencies of the environment modes (for 
the lowest non-trivial order in the strength of the spin- 
environment coupling) . The contribution of low frequen- 
cies <C B to <5b is exactly the field of a quadrupole, 
while high-frequency modes have the angular dependence 
of a quadrupole but a different radial dependence. 

Let us emphasize the novel points of our analysis in 
comparison to earlier work on BP in systems coupled to 
quantum or classical environments @> IS El El Some of 
these works have focused on the visibility of BP in spite 
of the dephasing, they did not find the modification of 
the BP. In it is absent because the Master equation 
used there neglects the effect of B ^ on the dissipative 
rates. A modification of the BP was found in , however 
the time-dependence of B(f) used there was too specific 
to see the geometric nature of this modification. Here we 
are able to consider a totally arbitrary (slowly varying) 
B(t) and show, for the first time, that the modification 
to the BP is geometric. We also observe that the effect 
of the BP on the dephasing rate is geometric. 

To find the BP we analyze the (directly observable) 
phase factor in the evolution operator for a given field 
dynamics B(t). In this phase we attribute to BP the 
contributions independent of ip [13 (those cx tp are part 
of the dynamical phase, and the terms 0(l/ip) are non- 
adiabatic corrections Jl3| ) . 

For open paths the BP is gauge-dependent 0], we 
show that this is not the case for the geometric dephas- 
ing, Im[<5$Bp] • The gauge-dependence of the BP is a 
consequence of ambiguity in the choice of instantaneous 
basis for a given B(t); B(t) defines the instantaneous 
z-axis but the i-axis may lie anywhere in the plane per- 
pendicular to B(t). Gauge transformations rotate the 
x-axis in this plane. The ambiguity is absent for closed 
loops, where the final basis must coincide with the initial 
one. In addition one can monitor the magnitude of the 
transverse spin component (dephasing); this magnitude 
is independent of the choice of x-axis, i.e. the dephasing 
(and corrections to it) is gauge-invariant even for open 
paths. 

Berry phase for an isolated spin. We evaluate the 




FIG. 2: The laboratory (x, y, z) and rotating (x, y, z) frames. 

BP for an isolated spin-half |l| using a rotating frame 
(RF)E1, the one we choose (see Fig. |2J has its z-axis 
along the field B(i) and eg _L e z (then ej(ip) = e^(0)). 
In this frame the magnetic field is B' = B + u>, where 
the angular velocity of the RF is 



P>e z 



(1) 



Transforming to the RF gives the Hamiltonian, H = 
UHspinU^ 1 +iUU~ 1 — — 7}(B + U))&, where the transfor- 
mation U = exp(i6a y /2) exp(ipa z /2) (with B, 8, ip be- 
ing spherical coordinates of B(i))[l||; B = (0,0, B) and 
ui = (u)x,u)y-,uJz) = ( — sin 9 6, cos 9 <£>) are coordinates 
in the RF. The evolution in the RF is a rotation of the 
spin about the z-axis by the angle § dt\R'(t)\. To lowest 
order in u> the accumulated phase is J dt{B(t) +u>g) [l7| . 
Subtracting the dynamical phase, we find the BP; 



<I>bp = / dtu- z {t) = <j) dip cos6> = <j)dBa(B) 



(2) 



In the chosen frame (gauge) the 'vector potential' a has 
only one non-zero spherical component a v — l/(_Btan#). 
Stokes' theorem is used to writing eq (0) as 



$bp= /rfSb(B), b(B) = V B x a(B) , 



(3) 



where Vb denotes derivatives w.r.t. the components of 
B, and the field b(B) is that of a charge-one monopole, 
b B = l/B 2 . 

Environmental contribution to BP. Adding a noisy 
environment-induced field X(t) = X(t)e z , the Hamilto- 
nian reads (in the lab frame) [7| 



H = -\B(t)a - \Xa z + H cnv (X) 



(4) 



This Hamiltonian models a situation in which the spin- 
environment coupling is strongly anisotropic. This is of- 
ten true for solid-state qubits, where various 'spin' com- 
ponents couple to entirely different environmental de- 
grees of freedom 0. [1 Our analysis is easily gen- 
eralized for multi-directional coupling^]. Below we 
express our results in terms of the statistical proper- 
ties (correlators) of the environment's noise, X(t). We 
consider any environment which gives rise to Marko- 
vian evolution, i.e. (X(t)X(0)) decays on a timescale 
t <C tp,Xdi ss = min(Ti,T 2 ,SB^ amh ) however we do not 
assume r <C B~ l . We further consider weak dissipation 



3 



H = -\B(t)a z - ±u}(t)a-±X(cos6a z -sm6a x )+H cnv . 

(5) 

We analyze this system in the eigen-basis of B{t)a z + 
u)(t)&, with the total field B'(t) w B(t) + uo z and the 
new angle 0' as defined in Fig. [21 

The small terms cx u>, (X 2 ) induce corrections to the 
evolution. As discussed above, the correction a w to 
the rate of phase accumulation yields the BP, while the 
correction oc (X 2 ) gives T-j - , , <5-E>Lamt>- Here we eval- 
uate the environment-induced BP cx (X 2 )lo. To this end, 
we reduced the problem with a time-dependent field to 
a problem with a stationary field (to leading order 
by going to the RF and then use the standard formalism 
for evaluation of dissipative effects. 

Consider the kinetic equation for the reduced density 
matrix of th e sp in p. Iteration of the Liouville equa- 
tion [HEl E3 HI El leads to a Dyson-type master 
equation: 

Pij(t) = \ [B'a^pi^+Jdtx^^it^Pi'Ah) • (6) 

The 'self-energy' £(i,ii) can be evaluated perturbatively 
in X\ diagrammatically it is the sum of irreducible di- 
agrams [22, E5| . For short-correlated noise r <C Th iss, 
one can make a Bloch-Redfield approximation [2(1 l2ll j 
giving Markovian evolution. We need only consider the 
off-diagonal density matrix element, pi2(t), as it contains 
all phases (and dcphasing) . After a secular approxima- 
tion (for r < B') Q, it is given by [H 

h 2 (t) = [iB'(t)+T(t)]p 12 (t), (7) 
r(f) ee fl^dt! £12,12(^1) exp B'(r)dr) .(8) 

We follow evaluating £12 12 to lowest order 

(2 nd -order) in X. This "golden-rule" calculation yields 



r(i) = - / dhSit-ti) cos0'(*)cos0'(*i) 

J —OO 

+^sme\t)sme'(t 1 )exp(-iJlB'(r)dr] , (9) 

where S(t - ir) ee (1/2) {{X (i)*(*i)> + (l(ti)X(t))} 
is the symmetrized environment correlator |26(. If we 
neglect all order-u effects, we obtain 

sin 2 9(t) 



r (t) 



2cos 2 0(i) 



O — -B + iO H 



fi + iO+ 



(10) 



where the integral is along the real axis and the function 
S(fl) is the Fourier transform of 5(r), it is real and even. 
The real part of Tq is the dephasing rate, which is given 
by the residue of the poles at O = 0, B. The imaginary 
part of To is the Lamb shift, coming from the principal 
value of the integral (only the sin 2 term contributes). 



Now we take into account the order-w term in B'(t) 
and the time-dependence of B'(£); variations of the angle 
between B and X and of B = |B| [i.e. 9{t 1 ) ^ d(t) and 
B(t) ^ B(t) in @]. To lowest order in u> the modifica- 
tion of the rate H1(J|) are: 



S v T(t) = 'up sin 2 9(t) cos9(t) x F(B) 
S e T(t) = i0 sin 9{t) cos 9(t) x G{B) , 
5 B T{t) = iB sin 2 9{t) x ^G'(B) , 



(11) 



the dot indicates (d/dt), the prime indicates (d/dB), and 



F(B) 
G(B) 



: S(0) 



dft S(O) (35 - 2fi) 



2^ 



2tt 5(ft 
5(0) 



B + i0+) 2 ' 
2S(fi) 



(12) 



(n-B + M+) 2 (o + io+) 2 



From Eq. Q) one sees that these C[w]-terms in T(t) 
generate O[w ]-terms in the total phase acquired by 
Pi2- $®bp = J[F{B) sin 2 9 cos 9 dip + G(B) sin 9 cos 9 d9 
+^G'(B) sin 2 9dB}. This is geometric, and corresponds 
to <5$ B p = / dB£a(B) (cf. eq. @) with the complex 
'vector potential' 



Sa v = B- 1 F(B) sin 9 cos 9, 
Sag = B~ 1 G{B) sin 9 cos , 5a B = ^G'(B) sin 2 



(13) 



For open-paths of B(t) this is our main result, the imag- 
inary part of this field (which gives geometric dephas- 
ing) is gauge- independent, while the real part is not. For 
closed paths of B(i) we use Stokes' theorem to write 
<5$ BP = JdS6b(B) (cf. eq. ©), finding that tfb(B) ee 
Vb x <5a has two non-zero components; 



Sb B = B- 2 F(B) (3 cos 2 0-1), 
Sbg = -B^F'iB) sin cos 0. 



(14) 
(15) 



Thus it is independent of Sag and 5a B , their only role 
(for closed paths of B(i)) is to form a 'pure gauge'; this 
is related to a symmetry discussed below. The angular 
dependence of Sh resembles that of a quadrupolc. For 
slow environment modes (fi -C B), Re[-F(-B)] oc B~ 2 and 
hence <5b(B) is a quadrupole field. For other environ- 
ment modes, ^b(B) has non-zero curl and zero diver- 
gence. Thus it is not a sum multipoles, it is the field 
generated by a pseudo-current, j(B) ee Vb x b(B), with 
one non-zero component 

j v {B) = B- 3 (6F{B) - B 2 F"{B)) sin0cos0. (16) 

Since we can ignore Sag and 5a b for closed paths, eq. i|13|) 
leads to the following pretty result: the total (complex) 
BP for a closed path is 



•Sbp = f dip 



— (B + SB^ + iTz 1 ). (17) 
di> z 



This result can be understood as follows. For a time- 
independent B the acquired phase is J dt£, where £ is 
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the term in parentheses in Eq. i|17|) . When B is time- 
dependent, B — > B + u> in the RF changing the phase by 
/ dt wV B £; with Eq. JT} this gives Eq. (|T7f) . 

Symmetry considerations. We now show that the 
vanishing of b v follows from the symmetries of the 
problem. Firstly, under time-reversal B(t),X(i) — > 
— B(— t), — X(— i) the BP between the excited and ground 
states is invariant, while in eq. © dS — » — <iS. This 
implies that b(— B, — X) = — b(B,Y). Secondly, if we 
instead rotate all spins and fields by angle n about the 
axis (p perpendicular to both B and X (e.g. B — > — B), 
the field b would also rotate. Thus, b(— B, —X) = 
i?Jb(B, Y). To satisfy both equalities, b v = 0. 

For isotropic coupling, the Hamiltonian's rotational 
symmetry guarantees that b v = bg = 0. The monopole 
(cf. eq. 10) is quantized so 5b b — 0; hence (5$bp = to 
all orders in the coupling to the environment. 

Contribution of the slow modes (O <C B). Here 
F(B) tx B~ 2 , so 5b is a quadrupole field. This effect 
of (for instance, classical) slow modes can be understood 
by noting that (B(t) + X(t)) varies adiabatically, hence 

$ B p - y"d(B+X)a m (B+X) = I dBa m (B+X) , (18) 

where a m is the field of a charge-one monopole at 
the origin. We disregard boundary corrections due to 
X(t P ) ^ X(0) and drop the term frfXa m (B + X) 
The remaining term is the field of a monopole at the 
point —X. In the multipole expansion, a m (B + X) = 
a m (B) + Y Q V Q a m (B) + iX ct X /3 V Q V /3 a ni (B) + •••. The 
first term produces the unperturbed BP, the second 
(dipole) term vanishes after averaging over the fluctu- 
ations [llj. Thus the quadrupole term gives the leading 
environment-induced modification of the BP. For noise 
along the z-axis this term reads fl a /3V a V^a m /6, with 
the quadrupole moment D a p = (X 2 ) diag(— 1, —1, 2). 

Concerning experiments. Our results imply that the 
traversal of the same loop B(t) in opposite directions 
would yield, apart from different BP's, different dephas- 
ing. We further note that a noise spectrum suppressed 
at = 0, B will increase Ti while having little effect on 
<5$bp, thus aiding the observation of <5$bp- 
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